CMB temperature trispectrum of cosmic strings 
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We provide an analytical expression for the trispectrum of the Cosmic Microwave Background 
(CMB) temperature anisotropics induced by cosmic strings. Our result is derived for the small 
angular scales under the assumption that the temperature anisotropy is induced by the Gott-Kaiser- 
Stebbins effect. The trispectrum is predicted to decay with a non-integer power-law exponent l~ p 
with 6 < p < 7, depending on the string microstructure, and thus on the string model. For Nambu- 
Goto strings, this exponent is related to the string mean square velocity and the loop distribution 
function. We then explore two classes of wavenumber configuration in Fourier space, the kite and 
trapezium quadrilaterals. The trispectrum can be of any sign and appears to be strongly enhanced 
for all squeezed quadrilaterals. 

PACS numbers: 98.80.Cq, 98.70.Vc 
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I. INTRODUCTION 



Although cosmic strings may be of various early uni- 
verse origins [H4To|. being line-like gravitational objects, 
they induce temperature discontinuities in the CMB 
through the Gott-Kaiser-Stebbins (GKS) effect [HG3. 
Direct searches for such discontinuities have been per- 
formed without success but do provide upper limits to 
the string tension U [l3|-[l5|. On the other hand, if cos- 
mic strings are added to the standard power-law ACDM 
model [16], it has been shown in Refs. [17|, ll8j that the 
CMB data are fitted even better if the fraction of the 
temperature power spectrum due to strings is about 10% 
(at I — 10). Such a fraction of string would even dom- 
inate the primary anisotropics of inflationary origin for 
I > 3000 With the advent of the arc-minute reso- 
lution CMB experiments and the soon incoming Planck 
satellite data, it is therefore crucial to develop reliable 
tests for strings [2(|, as to understand the non-Gaussian 
signals. The probability distribution of the fluctuations 
due to the GKS effect is known to be skewed and has 
a less steep decay than Gaussian a feature which 
can be explained in a simple model of kinked string [2l| . 
In Ref. [22J, we have studied the temperature bispec- 
trum induced by cosmic strings both analytically and nu- 
merically by using Nambu-Goto string simulations. We 
found good agreement between the analytical and nu- 
merical bispectrum for both the overall amplitude and 
the geometrical factor associated with various triangle 
configurations of the wavevectors. This agreement sug- 
gests that our analytical assumptions are capturing the 
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relevant non- Gaussian features of a string network and 
could be used to derive other statistical properties. In 
this paper, we present new results concerning the trispec- 
trum, i.e. the four-point function of the temperature 
anisotropy 23]. As pointed out in Ref. (22[, the bis- 
pectrum is generated only when the background space- 
time breaks the time reversal symmetry. Because our 
universe is expanding, the time reversal symmetry is in- 
deed broken and we get a non-vanishing string bispec- 
trum. On the other hand, the trispectrum can be gener- 
ated even in Minkowski spacetime and one may naively 
expect a stronger non-Gaussian signal than for the bis- 
pectrum. Motivated by this observation, we provide in 
this paper an analytical derivation of the string trispec- 
trum and study its dependency for various quadrilateral 
configurations in Fourier space. Given the fact that an- 
alytical predictions and numerical results exhibit good 
agreement both for the power spectrum [l^, [24| and the 
bispectrum [22j , we expect our result here to agree as well 
with the numerics. Performing such a comparison would 
however require a significant amount of computing re- 
sources which motivated us to leave it for a future work. 
Our main result can be summarized by Eq. (|54[) . Inter- 
estingly, the power-law behaviour of the trispectrum ex- 
hibits a non-integer exponent which can be related to the 
small scale behaviour of the string tangent vector correla- 
tor. In the framework of Nambu-Goto strings, this expo- 
nent is related to the mean square string velocity |25l. |26| 
and to the scaling loop distribution function [27ll28l]. The 
paper is organised as follows. In the next section we 
briefly recall the assumptions at the basis of our ana- 
lytical approach and derive the trispectrum in Sec. IIIII 
As an illustration, we apply our result to some specific 
quadrilaterals in Sec . 11 VI and exhibit some configurations 
that leads to a divergent trispectrum. They should pro- 
vide the cleanest way to look for a non-Gaussian string 
signal. 
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II. TEMPERATURE ANISOTROPY FROM 
COSMIC STRINGS 

In this section, we briefly review the general basics 
needed to calculate iV-point function of 9 = AT/T CMB at 
small angular scales. To study the correlation functions 
in the small angle limit, it is enough to consider 6 on the 
small patch of the sky. Then we can approximate this 
patch as two-dimensional Euclidean space, which simpli- 
fies the calculations. In this limit, the integrated Sachs- 
Wolfe effect generated by cosmic strings yields the tem- 
perature anisotropy in the light-cone gauge (24j 



k 2 O k = iek A J daX A {a)e lk - x{c 



where we have defined 



8tvGU, 



(1) 



(2) 



and X A [A = 1,2) is the two-dimensional string position 
vector perpendicular to the line of sight. We implicitly 
assume a summation on the repeated indices. It is now 
clear that the power spectrum, bispectrum, and higher 
order correlators can be evaluated in terms of correlation 
functions of the string network, as projected onto our 
backward light-cone. In order to evaluate the statistical 
quantities constructed over 0*-, the correlation functions 
of X A and X B have to be known. However, because 
0fc depends on X A and X B in a non-trivial manner, it 
is extremely difficult to derive meaningful consequences 
for the correlation functions without imposing additional 
conditions on the string correlators. In this paper, as 
done in Refs. [22|, [24j we therefore assume that both X A 
and X B obey Gaussian statistics, and this drastically 
simplifies our calculations. All the correlation functions 
of 0fc can now be written in terms of the two-point func- 
tions only. Using the same notation as in Ref. [2J|, the 
two-point functions of the string correlators are 



X A (a)X B (a>)) = ±6 AB V(a-a>), 
X A (a)X B (a')) = ^5 AB M(a-a'), (3) 
X A (a)X B {a')) = ^S AB T(a - a'). 



Note that an appearance of a term like e N(a — a') 
in the mixed correlator (X A X B ), where e AB is the anti- 
symmetric tensor with e 12 = 1, is forbidden due to the 
symmetry. As for the bispectrum, we also introduce the 
correlator 12211 



'yi 2 



T(a-a') = ([X(a)-X(a')] 



dci / da 2 T(ai - er 2 ) 



(4) 
(5) 



The leading terms are given by [2 
V(a)- 
I»- 

where we have defined 

£ = r'(oo), v 2 



v 2 a -> 
cr -> oo ' 

Per 2 a -> 
£er cr — > oo 



(6) 



(7) 



The correlation length £ is the projected correlation 
length on the backward light-cone, P is the mean square 
projected tangent vector (of order unity), v 2 is the mean 
square projected velocity (again of order unity). 



III. TEMPERATURE TRISPECTRUM 

In the flat sky approximation the four-point tempera- 
ture correlation function is defined as 



<e fcl e fc2 e fc3 e fc4 ) = T{k u k 2) fe 3 , fc 4 )(27r) 2 

x S(ki + k 2 + fe 3 + ki). 



(8) 



Using Eq. |T]) and a formal area factor A — (2ir) 2 5(0), 
the trispectrun^can be written as 



4 1 kik B k^k^ 
T(ki, fe 2 , fe 3 , fe 4 ) = e -t5aaSbb&cc5dd tjtjtjts 

>A yB t>C YDiS ab k a -Xt ' 



dcrida 2 dcr 3 dcr4 (XfXfXZXZe 

(9) 

with X A = X A (a a ), (a,b) £ {1,2,3,4} and fe a + fc 2 + 
^3 + ^4 = 0. We now assume Gaussian statistics and 
define 



pABCD yA y B t>C vD 
L = A t JL 2 A 3 A 4 7 

V = S ab k n ■ X h . 



(10) 

(11) 



The ensemble average in Eq. ^ can be expressed in 
terms of the two-point functions only 



(C ABCD e w ) = (C ABCD ) + (X A X B }(X^V)+ 



ABCD 
irA^x i <rB 



- A X B )(X^ 



(X A V)(X?V)(X^V)(X»V) 



1 

(12) 

where O denotes permutations of the labels {1,2,3}. 
Expressing X(a) in terms of X(cr) makes clear that 
all terms but the first involve the mixed correlators 



1 Notice that our denomination "trispectrum" here stands for the 
four-point function and contains the unconnected part. This one 
is however non-vanishing only for parallelogram configurations 
of the wave vectors. 
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(X(a)X(a')) . Since they give a small contribution com- 
pared to the others, we will not consider these extra terms 
in the following. A more detailed calculation would show 
that they induce corrections to the trispectrum scaling 
relative to the first term as 1/k and 1/fc 2 , respectively, 
and therefore are negligible at small angular scales. The 
trispectrum can be approximated as 



T(fei,fe 2 ,fe 3 ,fe4) - £ 4 ^S AA d BB S c5 S DB - 1 " ~ :i ' ; 



1L2 1,2 p 1,2 
1 2 3 4 



( (13) 

In terms of the two-point functions introduced in Sec.lHl 

(C ABCD ) = \s AB 5 CD V(a 12 )V(a 3i ) 

+ ls AC 6 DB V(a 31 )(a i2 ) + \s AD 5 BC V{a u )V{a 25 ) 1 

(14) 

where = a a — a^. As for (I? 2 ), replacing k^ with 
— ki — k 2 ~ k 3 , one gets 

(V 2 ) = ((fci • X u + k 2 ■ X 2i + k 3 • X 34 ) 2 ) , (15) 

where X a t, = X a — X^. As in Ref. [22j, one can show 
that 

(X u ■ X 24 ) = 1 [T(tru) + r(a 24 ) - r(ffia)] , (16) 

which can be used to transform Eq. (|15[) into a manifestly 
symmetric expression 



with 



K at> — k a ■ kj, 



(17) 



(18) 



From Eqs. (fl~3|) and (|14l) . we find that the integrand does 
not depend on 174 and the integration yields a factor equal 
to the total length of the strings L in the area A. In or- 
der to perform the integration over the other variables, 
one can again use the small angle approximation where 
all the k a are taken to be sufficiently large. The dom- 
inant parts then come from the small a length scales, 
the contributions from other regions being exponentially 
suppressed. This suggests we should Taylor-expand the 
two-point functions around a = 0. At leading order, us- 
ing Eq. ©, one gets fiy ~ Pa i4 a^ implying that (D 2 ) 
is a quadratic form in the variables 0^4. However, it ex- 
hibits a vanishing eigenvalue and the Gaussian integral 
cannot be extended to infinity since there is one direction 
of integration along which the exponent remains 
null. Let us notice that the situation is different than 
for the variable 04; the correlators are indeed function of 
such a flat direction whereas they do not depend on <t 4 . 
In order to get a sensible result, we therefore need to in- 
clude higher order corrections to the two-point functions. 

The behavior of T(a) at small scales is not trivial and 
many anal ytic al works have been devoted to its determi- 
nation I25LI25 129ll30j | . In the Polchinski and Rocha model 
of Ref. [25| , the next-to-leading order terms of the corre- 
lators (X ■ X) and (X ■ X) have a non-integer exponent. 
These correlators match with Abelian string simulations 
performed in Ref. [31] and can also be used to analytically 
derive the cosmic string loops distribution expected in an 
expanding universe. As shown in Ref. [28j], these results 
also match with the scaling loop distribution observed 
in the Nambu-Goto numerical simulations of Ref. [27| . 
As a result, we assume in the following a non-analytical 
behaviour for T(a) at small scales 



T(a) 



t 2 



2y 



(21) 



At this point, plugging this expression into Eq. ((TBI and 
performing the integrations along the lines done for the 
bispectrum is not possible (see Ref. [22]). Indeed, since 
the k a are forming a quadrilateral, contrary to the bis- 
pectrum triangle configurations, all the n a b cannot be 
positive thereby preventing some of the Gaussian inte- 
grals to be performed. 

We can nevertheless perform one integration by switch- 
ing to the more convenient integration variables 014, 024, 
(T34 and cr 4 . The Jacobian is unity and Eq. (|17j) can be 
rewritten in a non-symmetric form depending only on 
three of the variables: 



3 3 



(19) 



i=i j=i 



where 



fUj = 2 I r (^4) + I> j4 ) - i> l4 - an)] . (20) 



Notice that we are working in the light-cone gauge and 
therefore leave c\ and x as undetermined parameters 
since they cannot be straightforwardly inferred from the 
numerics performed in the temporal gauge. Nevertheless, 
because the correlation should be smaller as a becomes 
larger, c\ must be positive. Let us also mention the re- 
cent work of Ref. [26| suggesting that at very small length 
scales the correlator should become again analytic (with 
X = 1/2), i.e. that Eq. ([2Tj) would hold only for a > a c . 
However, as discussed in this reference, a c is shrinking 
with time in an expanding universe and at the times of 
observational interest, Eq. (f2"Tj) is expected to be valid on 
all the length scales we are interested in. 

With this next-to- leading order form of T(a), one ob- 
tains 



r(<r) ~ Pa 2 



Cl 



(i+ x )(i + 2x)e x 



(22) 
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and Eq. (IT9|) reads 



with 



/T) 2\ 1 ?2^i„ „ , „ K^cj)(a l4 ,a j4 ) 



"2(2 X + l)(2 X + 2)C 2 x 



where 



(23) 



<j>(a u ,a jA ) = (|a l4 | 2x+2 + M 2x+2 - |<r i4 - ^ 4 | 2x+2 ) . 

We can perform a linear coordinate transformation by in- 
troducing the set of orthonormal unit vectors (ei, e 2 , ei A 
and define three new coordinates xii X2 and \3 along 
these directions: 



(25) 



Xi = <5 u ( e i ■ fci)oj4, X2 = S tj (e 2 ■ ki)a j4 , 
X3 = e lJ '(ei ■ fci)(e 2 • fcj)<r w . 
We then have 

W = ^ + ,i) + 2(2x+1) g x+2)fr (26) 

X K J ^[cr l4 (x),cr i4 (x)] • 

The third coordinate %3 appears in (T> 2 ) only when the 
next-to- leading order terms in T(er) are taken into ac- 
count, which is consistent with the observation that there 
is a flat direction at leading order. The last term in the 
previous equation contributes little to the integrations 
over Xi and Xi- Hence we can safely say that the only 
non- vanishing component of x m the last term is X3- This 
is equivalent to include the next-to-leading order correc- 
tions only along the flat direction, i.e. for 

1 



°14 = J^fa ' fc ')( e 2 ' kj)X3, 



(27) 



where J is the Jacobian of the transformation given by 
Eq. (|25[) . Then, introducing the outer product coordi- 
nates by 



Wij = (ei • ki)(e 2 ■ kj) - (e x ■ k 3 )(e 2 ■ fcj) 



— ^ z \/^i^'j K ij i 



(28) 



2X+2 



one can show that 

2 

j2^+2 

1 

+ Ki| 2x+2 ) xl x+ \ 

(29) 

and other permutations. Finally, making use of identities 
such as 



«11^14,%) = -^T2 k >23\ 2X+2 xt X+2 , 



Kl2^(cri4,cr 2 4) = ^Tn ^12 (|W23| 2X+2 ~ |w 3 4 



«12 + «13 + «14 = 



one gets 



1 



Cl 



2 X +2 



(30) 



(31) 



y = -ki 2 \w 34: 



2x+2_ 







(32) 



Notice that Y 2 > for any quadrilateral because of 
the inequality (D 2 ) > 0. With these new variables, the 
trispectrum reads 



T(fci,fc 2 ,fc 3 ,fe 4 ) 



4Ak 2 k 2 k 2 k 2 J dXldX2 X J (33) 
{K 12 K 3i V[cr 12 (x)}V{a 3i {x)}+ O }e-2^ 2 >. 

At this stage, the Gaussian integrations over xi and X2 
are always finite, and for large enough wave numbers, 
i.e. kt; 3> 1, we can safely extend the integration range 
to infinity and also put xi = X2 = in V(aij). 
Eq. (j3"Tj) , the integration over xi and X2 yields 



From 



T(fci,fc 2 ,fe 3 ,fc 4 ) 



L 



+2 /Ipppp 



{K12K34W W34— \ V ( von— ) + O } 

2 X +2" 



(34) 



x exp 



with 



-c 2 Y* 



C2 



X3_y 

3) 



Cl 



2^x(2 x +l)(2 X + 2) 



(35) 



The integration over X3 ma y, a-priori, be performed in 
the same way. However, from Eq. (1321) . one can show 
that there is some particular configurations for which Y 
vanishes (parallelograms). As a result, one cannot push 
the integration up to infinity for those and one has to 
integrate only over the total string length. Notice that 
the integral depends on L only for the particular paral- 
lelogram configurations. As soon as Y 2 ^ 0, the small 
angle limit implies that Y 2 is large and the exponential 
function takes non- vanishing values only around vanish- 
ing X3- For this reason, we separate our analysis in two 
cases and first focus on the parallelogram case. 



A. Parallelogram configurations Y 2 — 

For parallelograms, the two opposite wavevectors form- 
ing the quadrilateral are anti-parallel and Y 2 strictly van- 
ishes. Without loss of generality, we assume k\ + fc 3 = 
and k 2 + fc 4 = 0. In this case, one has Wi 3 = w 24 = 
and we define 

w = wi 2 = W23 = k\k 2 sin (9 = fc 2 fc 3 sin#. (36) 

The integral in Eq. (|3"4")l can be evaluated along the flat 
direction X3/3, which is given by Eq. (|2"7) . 

CT14 = CT34 = W-y, cr 2 4 = 0. (37) 
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The integration range on Xil J is thus [—A, A] where 

L 



A 



2M 



(38) 



From Eq. (|34l) . the trispectrum simplifies to 
Tre 4 v 4 L 2 



T Q (k 1 ,k 2 ,k 3 ,k 4 



t 2 Akik 2 k 3 k 4 \w\ 
1 + 2COS2(0) ^/ 



i 2 Aklkl\siu9\' 



L/2 



V 2 (a) 



da 



(39) 



where we have neglected the integral in the last line. 
Since the correlator V 2 (a) is expected to be small at dis- 
tances larger than the typical correlation length £, this 
integral can be approximated by 



1 

1/2 



V 2 (a) 



da = -^V 2 (a )<2!«l, (40) 



Lv 



L 



where we have used the mean value theorem with a a G 
[0,£]. Under the scaling transformation k a — > bk a , the 
parallelogram trispectrum in Eq. (139[) scales as 



T (bk 1 ,bk 2 , bk 3 , bk 4 ) = b- 6 T Q (k 1 ,k 2 , k 3 , k 4 



(41) 



For parallelograms, it is important to recall that the 
trispectrum always gets a contribution from the uncon- 
nected part of the four-point function, which is purely 
given by a Gaussian distribution, 

T uc (fci, fe 2 , fc 3 , fe 4 ) - AP(k 1 )P(k 2 )+ • (42) 

As shown in Ref. [24] , the power spectrum is given by 

Lv 2 



P{k) = v 7 ^ 



Atk* 



(43) 



and the unconnected part of the trispectrum also scales as 
b~ 6 . Therefore the non-Gaussian contributions for paral- 
lelogram configurations remain of the same order of mag- 
nitude as the Gaussian ones, with the exception of the 
squeezed limit 9 — > 0. As we will see in the following, 
all other quadrilateral configurations have a scaling law 
which is different than Eq. (|41l) . 



B. Quadrilateral configurations with Y 2 3> 1 

In this case, the integrand in Eq. (|3"4")l takes non vanish- 
ing values only around \3 — an d we can safely extend 
the integration range over Xil^J to infinity, as we have 
done for \i an d Xi- One gets 



T oc (ki,k 2 ,k 3 ,k4 : ) ~ e 

x g(k 1 ,k 2 ,k 3 ,k 4 



-l/(2 X +2) 



fix) 

(44) 



The function f(x) is a number depending only on the 
power-law exponent x 



f(x) = 



-r 



i 



[4(2 X +l)( X + l)] 1/(2x+2) , 
(45) 



X + l V2X + 2 
and g({k a }) is the trispectrum geometrical factor 



g{k 1 ,k 2 ,k 3 ,k 4 



1 2 3 4 



-K12 [k\k 



2;,2 .2 \X+1 







-l/(2 X +2) 



(46) 



From Eq. ([44]) , we can derive various consequences worth 
mentioning. First, the sign of the trispectrum is com- 
pletely determined by the geometrical factor in Eq. (|46l) . 
which is manifestly symmetric under the permutation of 
two different wavevectors. The factor 



K12 ik\k\ 







(47) 



is always positive or zero. Also f(x) (f 0T a physically 
reasonable range of x) and c 4 are positive. Therefore the 
sign of the trispectrum is given by the factor 



K12K34 + K13K24 + K14K23, 



(48) 



which can be positive or negative according to the quadri- 
lateral under scrutiny. 

Secondly, under the scaling transformation k a — > bk a , 
the geometric factor scales as 



g(bki,bk 2l bk 3 , bk 4 ) = b p g{k 4 ,k 2 , k 3 ,k 4 ), 



with 



P = 6 



X + l 



(49) 



(50) 



Contrary to the case of the power spectrum, of the bis- 
pectrum, and of the parallelogram configurations, the 
leading term of the trispectrum scales with a non-integer 
power-law exponent. For x > 0, the trispectrum decays 
slightly faster at small scales than the bispectrum. Let 
us recap that in the temporal gauge, the string tangent 
vector correlation function exponent x is a small quan- 
tity related to the expansion rate of the scale factor and 
to the mean square velocity of strings (25[. This is cer- 
tainly also the case in the light-cone gauge and one may 
be able to use the trispectrum to distinguish between 
different models of strings. For instance, in Abelian 
Higgs numerical simulations, the strong back-reaction in- 
duced by scalar and gauge radiation produces a mean 
square velocity lower than in classical Nambu-Goto sim- 
ulations [3J [Hj]. Meanwhile, the loop distribution ob- 
served in Nambu-Goto simulations has a power-law ex- 
ponent which is uniquely given by x [271 ] . Interestingly, 
the scaling exponent is different from the one associated 
with parallelogram configurations. These two different 
scaling behaviors may actually be used to distinguish the 
trispectrum by cosmic string with the one generated by 
other sources. 
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C. Interpolating trispectrum for all quadrilaterals 

When Y 2 ~ but non- vanishing, i.e. for quadrilat- 
erals close to parallelograms, one cannot push the inte- 
gration range in Eq. (|34[) to infinity. Contrary to the 
case Y 2 = 0, the integration over x?>l J cannot be per- 
formed explicitly in this case. Nevertheless, we can make 
some approximations. First, for configurations close to 
parallelograms, two of the quantities are expected 
to be small, say w\ 3 and u>24- For those, one can re- 
place the V(a) functions in Eq. (|34|) by v 2 . On the other 
hand, one expects the other w m n factors to be large and 
Eq. ([3"^]) has terms involving the product V(a)V(a'). As 
for the parallelograms, we expect those to be at most 
of the order v 4 £/L, which can be neglected compared to 
the terms in v A . With an integration range over x 3 jj 
given by [—A, A], where A(fe Q , L) has still to be specified, 
performing the last integration over \3 yields 



4 v L£ / \ 



-l/(2 X +2) 



fix) 



1 2 3 4 



where 7 n (a, x) denotes the normalised incomplete lower 
gamma function 



7„(a,x) 



~ r(a) 



(52) 



In the limit Y 2 0, this expression matches with 
Eq. <[39]) for A = L/(2\wi 2 \). In order to interpolate 
between Eqs. (|39p and (|44l) we can replace the geomet- 
rical factor in Eq. (|5ip by the factor g({k a }) and chose 
the cutoff A to be 



A 



2L 



1 1012 1 + |tWl3| + |«>14| + \w 2 3\ + \w2i\ + \tVu\ 



(53) 



K12K34 + K13K24 + K14K23 

Our interpolation formula for the trispectrum finally 
reads 



T(fc 1 ,fc 2 ,fe 3 ,fc 4 )^e 4 ^(c 1 C 2 ) 



-l/(2 X +2) 



fix) 



x 7„ 



1 



2X + 2' 



,c 2 Y 2 A 2 ^ +2 )g(k 1 ,k 2 ,k 3 ,k 4 ), 



(54) 



with A given by Eq. ([S3]) and g{{k a }) by Eq. (gBJ). For Y 2 
large, the gamma function is close to one and we recover 
Eq. (|4~4"]) . The limit Y 2 = gives again the leading order 
of Eq. J3JJ. 



IV. GEOMETRICAL FACTORS OF 
SYMMETRIC QUADRILATERALS 

In this section, we explore the dependency of the 
trispectrum geometrical factor given by Eq. (|46p for some 





FIG. 1: Quadrilateral configurations for the trispectrum 
wavenumbers. The left is referred to as the "kite" quadri- 
laterals with two adjacent sides of equal length and the other 
two sides also of equal length. The right quadrilateral is a 
"trapezium" which is defined to have two opposite sides par- 
allel. 



symmetric quadrilateral configurations of the wavevec- 
tors. 



A. Kite configurations 

Let us first consider a quadrilateral like the one given 
in Fig. □ (left). From Eq. ([32]), one gets 



Y 2 =k 6+ ^y 2 (d,a) 



(55) 



with 



i+x 



y 2 (6,a) = [sin 2 (6/2)] 



x { 2sin(0/2) — i x 



cos a 



cos 2 (a- 6»/2) 



i+x 



- 2sin(0/2) 



sin(a + 0/2) 



cos a 



cos 2 (a + 0/2) 



+ 4 1+x sin 2 (0/2) [cos 2 (0/2)] 



cos z a 
i+x cos(2a) 



i+x 



cos 2 (a) 



4 1+x cos(0) [sin 2 (0/2) tan 2 (a)] 



i+x j 

From Eq. ([46]) . the geometrical factor reads 

cos 2 (a) [1 — 2 cos(2a) cos( 



g(k 1 ,k 2 ,k 3 ,k 4: ) 



sin 2 (0/2) 



1 



(56) 



(57) 



J t p y 2/(2+2 X ) 



As expected from the trispectrum scaling law, the kite 
trispectrum decays as l/k p at small angular scales. The 
overall amplitude is however amplified for squeezed con- 
figurations and diverges for 0—^0. For small, the 



7 



30 
25 
20 
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FIG. 2: Trispectrum geometrical factor for the kite quadrilat- 
erals as a function of 0, plotted for various values of a. The 
trispectrum is enhanced in the squeezed limit 9 — > 0. The 
bump for 6 P = 7r — 2a corresponds to the parallelogram limit 
for which the unconnected part is no longer vanishing. 



leading terms of the previous expression are 



8 cos 2 (a), „ „ . 
g , - V l-2cos2a 



x J 2(1 + x) tan 2 (a) - 1 + 4 X (1 - tan 2 a) 



-l/(2x+2) 



(58) 

The sign of the kite trispectrum is the same as 1 - 
2 cos(2a) cos(0) and, at small 0, is negative for a < tt/6 
and positive otherwise. As for the bispectrum, we recover 
that squeezed configurations are the most sensitive to a 
string signal, certainly due to the elongated temperature 
discontinuities induced by the GKS effect. In Fig. [21 we 
have represented the full geometrical dependency coming 
from Eq. (l54l) as a function of and for various values 
of a. For convenience, we have chosen \ — 0.29, c 2 = 1, 
k = 1 and L = 20. The incomplete gamma function con- 
tributes for configurations close to the parallelogram ones 
which appear as a bump in Fig. [2] for p = tt — 2a. For 
the kites, the argument of the gamma function simplifies 
to 



c 2 F 2 A 2 x+ 2 = k 2 



2L 

T 



2X+2 



y 2 (9,a) 



2(2 X + l)(2 X + 2) 
[l-2cos(2a)cos(0)]~ 2(x+1) 

2(X+1) 



{2 sin(0) + [cos(0) - 1] tan(a)} 

(59) 

As can be seen on this plot, we recover the change of 
sign when a crosses the value ir/6. The bump at # p = 
tt — 2a corresponds to the parallelogram limit of the kite 
configuration for which y 2 (9, a) — > 0. 



B. Trapezium 

Let us next consider a quadrilateral given by the right 
side of Fig Q] having two opposite sides parallel. Without 
lost of generality, one can assume that the upper side is 
of smaller length than the bottom. Denoting their ratio 
by sin 2 (/3), after some algebra, the factor Y 2 is still given 
by Eq. (|5H} with 



y 2 (9,a) 



[sin 2 (0)] 
s 2 (/?)] 2x+1 



x+i 



sin 2 (a + 6 
sin 2 (a) 

sin 2 (/3)] 2X+1 



X+2 



(60) 



l2x+2 



tan 2 (/3) [sin 2 (/?)]" 
Similarly, the geometrical factor reads 

sin(a) sin(0) — 3 cos(a) cos(0) 



g(k 1 ,k 2 ,k 3 ,k 4: ) 



sin(a) 
sin 2 (6») 



sin 2 (a) 



-I (p-3)/2 



sin (a + i 



{l-[cos 2 (/3)] 2x+1 -[sin 2 (/3)] 2x+1 } 



kP 

in 4 (/3) [tan 2 (/3)] 

l/(2x+2) 



(p-4)/2 



(61) 

As expected, the trapezium trispectrum decays with the 
power-law exponent k~ p . The overall amplitude is again 
amplified for elongated configurations and diverges for 
9^0. For convex quadrilaterals, assuming < 9 < 
■k — a, the sign of the trispectrum is given by the first 
term of Eq. (|6"Tj) . As a result, it is is negative for 9 < 9 S 
and positive otherwise, where 9 S is given by 



sin(a) 



y^9cos 2 (a) + sin 2 (a) 



(62) 



For an isosceles trapezium with a = 9, the change of sign 
occurs at 9 S = tt/3. Finally, in Fig. [31 we have plotted 
the full geometrical dependence as a function of 9, for 
various values of a. For the trapeziums, the argument of 
the gamma function is 



C2 y 2 A 2 *+ 2 = fc 2 



2L\ 



2X+2 



2(2 X + l)(2 X + 2) \ i J 
x [sin(a) sin(0) - 3cos(a) cos(6»)]~ 2(x+1) 



y 2 (9,a) 



sin(0) sin(0 + a) [3 • 



-2(x+l) 



sin(a) sin (j3) 



(63) 

The divergence for the parallelograms visible at = tt — a 
comes again from the squeezed shape. Imposing a fixed 
value of sin 2 (/3) implies that such parallelograms are in- 
finitely elongated. The configuration with > tt — a 
are self-intersecting trapeziums having a butterfly shape. 
Their squeezed limit occurs for — > tt for which the 
trispectrum is again strongly enhanced. 
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FIG. 3: Trispectrum geometrical factor for the trapezium 
quadrilaterals as a function of 9, plotted for various values of 
a. For convenience, the ratio of the two parallel sides has been 
fixed to 3/4 and \ = 0.29. The divergence in the squeezed 
limit occurs at 9 — ► but also at 9 P — n — a for infinitely 
elongated parallelograms. For 9 > 9 P , the trapeziums are no 
longer convex and represent "butterfly" configurations which 
are squeezed for 9 — > n. 



V. CONCLUSION 

In this paper, we have analytically derived the CMB 
temperature trispectrum induced by cosmic strings us- 
ing the string correlation functions in the Gaussian ap- 
proximation. The trispectrum generically decays with a 
non-integer power-law behaviour at small angular scales 
which depends on the string microstructure through the 
behaviour of the tangent vector correlator on small dis- 
tances. Its eventual detection and measurement may 
therefore help to distinguish between different string 
models. We have also found that the trispectrum di- 
verges, in the framework of our approximations, on 
all squeezed configurations whose measurements remain 
however limited by the finite experimental resolution. In 



fact, such a non-integer power-law is linked to the exis- 
tence of a "flat direction" at leading order and the four- 
point function ends up being sensitive to the next-to- 
leading order string tangent vector correlator. This sit- 
uation is also present in the n-point function and we do 
expect all of the higher n-point function to exhibit non- 
integer power-law behaviours. Since this situation was 
not encountered for the two- and three-point functions, 
the next step will be to compare our results here with 
the trispectrum computed from CMB maps obtained by 
string network simulations. 

Finally, let us notice that we have not attempted to 
make any comparison with a CMB trispectrum produced 
by primordial non-Gaussianities of inflationary origin. 
The situation is nearly the same as it is for the string 
bispectrum (22[. The so-called tnl arid gNL parame- 
ters quantify the amplitude of the primordial four-point 
function of the curvature perturbation on super-Hubblc 
scales. As a result, the induced trispectrum of the CMB 
temperature fluctuations strongly depends on the CMB 
transfer functions and exhibits damped oscillations with 
respect to the multipole moments. Here, we have dire- 
cly derived the CMB temperature trispectrum produced 
by the strings and it would therefore make no sense to 
find an associated tnl and <?nl ■ An alternative approach 
might be to estimate what values tnl and c/nl would 
assume in a primordial-type oriented data analysis if 
the non-Gaussianities were actually due to strings. This 
could be done with a Fisher matrix analysis for a given 
experiment but we leave this question for a forthcoming 
work. 
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